Multi-pass classical vs. quantum strategies in lossy phase estimation 
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The use of classical multiple-pass approach for phase estimation which mimics the behavior of 
the NOON states, is compared with quantum techniques. It is shown that in the presence of losses 
its performance is significantly worse than the one of the optimal quantum strategy. 
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I. INTRODUCTION 

Optical interferometry [1] is capable of ultra pre- 
cise measurement of various physical quantities such as 
length, time, temperature, velocity provided they can be 
translated into a phase shift ip of the light beam. In an 
interferometric experiment, light intensities detected at 
the output ports carry information on the true value of 
ip. Once intensities are measured the estimation proce- 
dure may be carried out resulting in an estimated value (p. 
The intrinsic character of the detection process, however, 
is stochastic. Within the semi-classical approach this is 
attributed to the quantum nature of atoms constituting 
the detectors, which under the interaction with classical 
light field are being excited. The Poissonian statistics 
governing the number of excited atoms, and therefore 
the number of detector clicks - the shot noise - leads to 
phase estimation uncertainty scaling: 5ip cx \j\ffi where 
n is the mean number of detector clicks. 

Quantum description of light, provided by quantum 
optics, allows for a deeper analysis of what are the ulti- 
mate limits of phase estimation precision. Semi-classical 
Sip cx l/\/n scaling can be recovered when one restricts 
oneself to the use of coherent states of light, and n is iden- 
tified with the mean number of photons used. What is 
more interesting, however, is that other inherently quan- 
tum states of light, such as e.g. squeezed states, twin- 
Fock, NOON states, offer in certain scenarios quadratic 
improvement of phase estimation precision leading to the 
Heisenberg limit (HL): 8ip cx 1/n [2-9]. The quantum- 
enhanced precision comes at a price. Quantum states 
achieving the HL require highly nontrivial preparation 
procedures, such as strong squeezing or large entan- 
glement. Moreover, the ultra-high phase sensitivity of 
the states is usually accompanied by low robustness to 
environmental decoherence. In particular, photon loss 
renders the NOON states useless and other more robust 
classes of quantum states needs to be employed [10-18]. 
Fragility of quantum states and difficulties in their prepa- 
ration are the challenges that need to be met if quantum 
enhanced interferometry is to be implemented in prac- 
tice. 



The situation prompted ideas on overcoming the shot 
noise limit without resorting to the use of sophisticated 
quantum states, but rather allowing some minor modifi- 
cations in interferometric experiment. One of the most 
simple techniques is the multi-pass approach, in which 
a light beam is allowed to pass the sample many times 
and thus acquire a multiple of the phase to be estimated 
- in this way it mimics the behavior of the NOON state 
[19, 20]. Similar ideas were also analyzed in the context 
of precise clock synchronization [21, 22]. Obviously such 
an approach can be seen as a change in the rules of the 
game, and the use of quantum states may be still argued 
to be the only way to reach the HL. 

The presence of losses, complicates the relation be- 
tween classically based and quantum based strategies 
even further, neither NOON nor multi-pass approach is 
the optimal strategy. In what follows the optimal use 
of classical states for interferometry in the presence of 
losses is investigated and compared with the optimal use 
of quantum states. It is shown that no straightforward 
classical multi-pass technique can reach the precision of- 
fered by quantum states when losses are present. 




FIG. 1: The generic Mach-Zennder interferometer fed with a 
classical light beam. The phase shift tp inside the interferome- 
ter modulates the output intensities. Possible losses accompa- 
nying the phase shift are represented by a beam splitter with 
power transmission r\. In multi-pass strategies light passes 
the sample k times. 



II. CLASSICAL INTERFEROMETRY IN THE 
PRESENCE OF LOSS 
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A paradigmatic model for the study of phase estima- 
tion is the Mach-Zehnder interferometer (MZI) depicted 
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in Fig. 1. A classical light beam is fed into the lower input 
port of the interferometer. The beam is divided at the in- 
put beam-splitter (transmissivity T), experiences a rela- 
tive phase delay tp caused by a phase shift element placed 
in the upper arm, and — if no multi-passes are allowed — 
is recombined at the output beam-splitter (transmissivity 
50%). Finally, its intensity is being measured at the two 
output ports. The phase shift element, apart from intro- 
ducing a phase delay, typically is also a source of losses. 
This can be modeled by placing an additional beam split- 
ter with transmissivity r\ in the upper arm. A change of tp 
can be sensed by observing the change of light intensities 
at the output ports. 

The probability that the two detectors register n\ and 
712 clicks, provided the unknown phase shift equals tp can 
be expressed as a product of two Poisson distributions 
p( n i, n 2\tp) = Pni ( n i)Pn 2 { n 2), with respective mean val- 
ues 

fi\ = An (1 — v cos tp) , t%2 — Ah (1 + v cos tp) (1) 

where h is the mean number of photons in the input 
beam, while: 



A 



i-T(i-n) 



2 v /r(i-?> 
i-T(i- v ) ' 



(2) 



For T] = 1 (lossless case), the choice T — 50% causes 
the functions h\{tp), h 2 {tp) to enjoy both the v = 100% 
visibility and the maximum amplitude A = h/2. Any 
other value of T will lead to a decreased visibility and 
hence it is no surprise that T = 50% is the optimal choice 
for an interferometric experiment. For n < 1, however, 
the optimal choice for T is no longer obvious and one 
needs more precise tools to find the solution. 

Knowledge of p{n%, n 2 \tp), and measurement results m, 
n 2 allows for an estimation of the value of tp. One of the 
key tools in estimation theory is the Cramer- Rao inequal- 
ity [23] which provides a fundamental bound for the root 
mean square error dtp of estimation expressed in terms of 
the Fisher information F: 



Sip > Stpn 
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dp(n 1 ,n 2 \tp) 



dtp 

(3) 

Loosely speaking, the larger is the derivative of 
p(ni,n 2 \tp) over tp, the larger is the Fisher information, 
the lower is the achievable variance of estimation and con- 
sequently the better is the estimation precision. For the 
setup considered the Fisher information reads explicitly: 



F 



4nAv 2 sin tp 
- v 2 (cos 2tp + 1) 



(4) 



It is intuitive to expect that in the presence of losses 
(?7 < 1), the optimal choice for T should correspond to 
more light being transmitted into the lossy arm in or- 
der to keep the visibility of interference high. Choos- 
ing T = 1/(1 + 77) we retain 100% visibility (v = 1), 




FIG. 2: Mean number of clicks registered by detectors (top) 
for r\ = 0.1 when the transmission of the input beam splitter of 
Mach-Zehnder interferometer is chosen to maximize visibility 
(dashed), or to minimize estimation uncertainty (solid). The 
dependence of uncertainty on tp (bottom) illustrates the fact 
that the maximum-visibility strategy works equally well for 
the whole 2-r region, whereas the optimal strategy requires 
operating in the vicinity of phases indicated by arrows, where 
it achieves the standard interferometric limit. 



and get F — 2ray/(l + n). Surprisingly, this is not the 
optimal choice. The optimal choice T = 1/(1 + ^/rj) 
results from a compromise between keeping visibility 
high and loosing as little light as possible, and yields: 
F — 4nrysin 2 tp/(l + f] — 2^Jr\ cos lip). The price to pay 
is the dependence of F on the value of tp, and the opti- 
mally is achieved only in the most sensitive estimation 
■ region corresponding to tp — ir/2,3ir/2, where hi(tp), 
h 2 (tp) curves are the steepest. In this region, we get 
F = Ahrj/ (1 + y/rj) 2 and the minimal uncertainty: 



Sip: 



SIL 
mill 



lnr\ 



(5) 



The above classical precision obtained for the optimal 
choice of beam splitter transmission is referred as the 
standard interferometric limit (SIL) [16, 17], with its 
characteristic \j\fn~ shot-noise scaling. Fig. 2 depicts 
the differences between the optimal and the maximum- 
visibility estimation strategy. 

The intuitive argument: the steeper the curves the bet- 
ter the precision, should be taken with care. In partic- 
ular for the 100% visibility case, Fisher information is 
constant over the whole [0, 27r] region even though the 
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steepness of the curves varies greatly. When the curves 
become flat, one would expect vanishing of the Fisher in- 
formation. This is not the case, however, since simulta- 
neously one of the mean values n\{ip), fi2(ip) approaches 
zero, and this causes the l/p(n\, ri2\ip) term to blow up 
keeping the Fisher information constant. 



III. QUANTUM STRATEGIES 

Optimal quantum strategies allow for a significant im- 
provement over the SIL even in the presence of losses 
[15-17]. For the sake of generality one replaces the in- 
put beam-splitter of MZI with a general quantum state 
preparation device, and analogously the output beam- 
splitter plus detectors with a general quantum measure- 
ment. In the absence of losses (77 = 1), the optimal n 
photon state that should be prepared by the preparation 
device is the NOON state [8] \ip) = (|n)®|0)+|0)<g)|n})/\/2, 
where the first (second) factor in the tensor product cor- 
responds to the number of photons in the upper (lower) 
arm of the interferometer. The minimal estimation un- 
certainty corresponds to the HL: 



Sip 



HL 
min 



1 

n 



(6) 



which is a quadratic improvement over SIL. The improve- 
ment stems from the fact that the two superposition 
terms acquire a relative phase factor exp(inip), which 
causes measured fringes to be n times denser (and hence 
n times steeper) than if a single photon or classical light 
were used. In the presence of losses, however, NOON 
state performance deteriorates rapidly. Even, if we con- 
sider a NOON state with optimally distributed weights 
among both superposition terms the increase of uncer- 
tainty is exponential in n: 



Sip 



NOON 
min 



1 



V 



i/2 



2 nr] n/2 



(7) 



and for large n its advantage over classical strategies 
is quickly washed out even in the presence of moder- 
ate losses. Exponential growth of uncertainty may be 
avoided if instead of a single n photon NOON state, a k 
photon NOON state is sent m — n/k times through the 
interferometer [16, 17] 
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Relaxing the integer constraint on m and fc, the opti- 
mization over the choice of k <G [l,n] yields 



Sip' 
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where the optimal choice of k for the three regimes indi- 
cated above is k = 1, k = 1.478/| I1177I (i.e. the solution 
of 1 + r/ fe / 2 + k In T] = 0) and k = n. 

More robust states have been proposed [16, 17] which 
outperform both classical and NOON state based strate- 
gies and raise hopes for beating SIL in realistic environ- 
ments. The states are found by numerical optimization 
of quantum Fisher information which unlike its classical 
counterpart in Eq. (3), does not require specifying the 
actual measurement performed on quantum states. For 
a general state of n photons distributed among the two 
arm of the interferometer 



^ a s \s) ® |n - s) 



(10) 



s=0 



phase estimation uncertainty reads: 




\ " 



\Z n s=iXssBf(T,)f 



-1/2 



(11) 

where x s = \a s \ 2 , Bf(jy) = (^"'(l - vY- Optimal 
states result from minimization of Sip^ in over x s . These 
states have complex structure and their optimality is the 
result of a subtle interplay between the need for quan- 
tum enhanced phase sensitivity and robustness to losses. 
While the estimation uncertainty is greatly reduced when 
using these states for realistic system parameters the HL 
cannot be retained [16, 17]. 



IV. MULTIPLE-PASS CLASSICAL STRATEGIES 

Comparing different estimation strategies requires a 
precise definition of resources that are available. In most 
application the mean number of photons used is a rel- 
evant quantity. Other approaches are conceivable how- 
ever. One could fix the mean photon number traveling 
through the sample, which is relevant in situations when 
not the total power consumed is the limitation but only 
the fraction of it that is absorbed by the sample. One 
of the examples in which this distinction is relevant is 
the multi-pass approach to phase estimation. Instead of 
sending light only once through the sample, one may re- 
flect it and make it pass the sample k times, acquiring 
a multiple kip of the estimated phase. This simple trick 
mimics the behavior of NOON states, making interference 
fringes k times denser and therefore allowing for better 
phase sensitivity. According to Eq. (3), Fisher informa- 
tion depends quadratically on the derivative of the prob- 
ability distribution over phase. Increasing the density of 
fringes k times makes the Fisher information grow as k 2 , 
and consequently results in 1/k scaling Sip m i n . 

Two separate cases may be considered: when multiple- 
passes are for free or when they are included in the overall 
budget as as a resource equivalent to photons. 
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A. Multiple-passes as a resource 

If the total resources are to be fixed to n, the intensity 
of light passing through the sample k times has to be 
reduced to n/k. This way we obtain interference fringes 
which are k times denser, at the cost of replacing n with 
n/k and n with n k . Mathematically, this strategy is 
equivalent to the NOON chopping strategy with k photon 
NOON states sent through the interferometer n/k times, 
resulting in Eqs. (8). In the context of clock synchroniza- 
tion a analogous equivalence in susceptibility to external 
decoherence has been observed between quantum GHZ 
states and multi-pass strategies using single qubit states 
[21, 22]. 



B. Multiple-passes for free 

If the number of bounces k does not contribute to the 
consumed resources, we can choose the value of k that 
yields the highest sensitivity of estimation. Mathemat- 
ically, this situation is equivalent to replacing ip — > kip, 



n k . Choosing the optimal T = 1/(1 + n k / 2 ) in he 



MZI we get the estimation uncertainty bounded by: 



<Vmin = 



nk/2 



(12) 



Notice that in the absence of losses the precision can 
be arbitrary good since the bigger the k the lower the 
uncertainty. When losses are present, minimal Sip min is 
achieved for k = —2(1 + £)/\nn, where £ w 0.278 is the 
solution to £e^ +1 = 1. The resulting uncertainty: 



x „mp 

<Vmin = 



In n 
4Vl£' 



(13) 



V. MULTIPLE-PASS QUANTUM STRATEGIES 

For the full picture we also need to investigate the use- 
fulness of multi-pass approach when employing quantum 
states. 



A. Multiple-passes as a resource 



e lks f . The same effect, however, can be obtained with 
a single pass of ks photons prepared in state: 



5> s 

s=0 



ikscp 



\ks) (g> \kn — ks). 



(15) 



Therefore, there is no advantage, except for technical dif- 
ficulties in preparing multi-photon quantum states, in us- 
ing multiple-pass approach if multiple-passes are treated 
as a resource equivalent to photons. The same conclusion 
holds also in the presence of losses, since the probability 
of loss behaves in exactly the same way, irrespectively 
whether a given number of photons pass k times through 
the sample or a k- multiple of them pass through the sam- 
ple once. 



B. Multiple- passes for free 

The above argument clearly cannot be applied in the 
case when we may freely choose the number of passes and 
do not account for them in the total amount of resources 
consumed. This approach will in general improve esti- 
mation precision as compared with the optimal quantum 
single-pass strategy. When the general n photon state 
given in Eq. (14) passes the sample k times, we can sim- 
ilarly as in the single pass case derive the formula for 
minimal uncertainty and the only difference in the final 
result is the replacement of n by n k and the division of 
the whole formula by k: 



Sip' 



QMP 
min 




-1/2 



(16) 

The optimal strategy can be found by direct minimiza- 
tion of Stp^n over x s and k parameters. Due to the 
additional parameter k which enters the formula in a 
nontrivial way, numerical optimization is a bit more com- 
plicated than in a single pass case. Unlike S(p^ in from 
Eq. (16) which was proven a convex function of x s [17], 
S'P^f^' is not a convex function of k. One needs to take 
an additional care not to end up in a local minimum in- 
stead of a global one. 



Neglecting losses for the moment, consider a general 
n photon quantum state which has experienced a phase 
shift <p in the upper arm with respect to the lower arm 
of the interferometer: 

n 

\i/>) =^a s e- is v|s)® \n-s). (14) 

Multiple-passes give the advantage of multiple acquisi- 
tion of phase and hence better sensitivity. If the light 
passes the sample k times, we need to replace e~ ls '- p with 



VI. COMPARISON 

Comparison of estimation uncertainties achievable us- 
ing all strategies presented in the paper for a particu- 
lar value of interferometer transmission n — 0.6 is de- 
picted in Fig. 3 as a function of number of photons used. 
Quantum single-pass strategy clearly outperforms clas- 
sical single-pass strategy as well as classical multi-pass 
strategy where multiple-passes are treated as a resource. 
Notice the \/\fn character of the curves for classical 
strategies, as opposed to a steeper one for the quantum. 
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FIG. 3: Log-log scale plot of phase estimation uncertainty for 
transmission n = 0.6 for optimal classical single-pass strategy 
(black, dotted), quantum NOON chopping strategy or equiv- 
alently classical multi-pass strategy (black, dashed), optimal 
quantum strategy (black, solid), where n denotes total re- 
sources consumed equal to the number of photons used times 
number of passes through the sample. If multiple passes are 
not considered a resource, the uncertainty decreases, which 
is depicted by two gray curves, representing optimal classical 
multi-pass strategy (gray, dashed) and quantum multi-pass 
strategy (gray, solid). Lower shaded area corresponds to un- 
certainties smaller than that given by the Heisenberg limit: 
1/n, while the upper shaded area corresponds to uncertain- 
ties larger than the shot noise limit 1/y/n. 

More importantly, with growing number n the advan- 
tage of the optimal quantum strategy is more evident. 
Nevertheless, due to losses it is not powerful enough to 
provide 1/n HL which is achievable with NOON states in 
the absence of losses. 

When multiple-passes are allowed for free, the uncer- 
tainties of both the optimal classical and quantum strate- 



gies are significantly decreased, falling for small n even 
below the HL. Notice however, that the slope of the clas- 
sical curve manifests again its 1 / ^fn character. Quantum 
curve is again steeper and proves its advantage over clas- 
sical one more and more the bigger n becomes. 



VII. SUMMARY 

This paper analyzed in detail the impact of multiple- 
passes for phase estimation using both classical and quan- 
tum strategies. In particular it has been shown that while 
in the absence of losses classical multi-pass strategy may 
in some sense be regarded as equivalent to the optimal 
quantum strategies this is no longer the case when losses 
are present. It might sound as a paradox, since we are 
used to the fact that quantum states manifest their ad- 
vantage over classical ones best in well controlled envi- 
ronment with low noise. In the phase estimation prob- 
lems, however, if multiple-pass approach is allowed then, 
indeed it is only in the presence of losses that quantum 
states retain their advantage over classical based strate- 
gies in the quest for precise phase estimation. 
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